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B meson mixing at CDF II 

M. Rescigno (for the CDF CoUaboration)'' 

'^Istituto Nazionale di Fisica Nucleare, Sezione di Roma, P.le Aide More 2, 00185 Roma, Italy 

We present the first limit on mixing frequency obtained using 360 pb^^ of Tevatron Run II data with the 
CDF II detector. We derive Am^ > 7.9 ps"^@ 95% C.L. with a sensitivity of 8.4ps"^ 
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1. Introduction 

In the Standard Model (SM) flavor mix- 
ing of neutral mesons is induced by second- 
order charged weak interactions through box di- 
agrams. For B° the mixing frequency, Ato^, is 
a /s^-BB^IVtdp and thus provides a constraint 
on the usual unitarity triangle whereas for i3° 
Arris C!C i^s, |Vtsp and is much higher since 
\Vts\ » \Vtd\- Theoretical uncertainty on the 
hadronic parameters make the constraint from 
the precise Arrr^ measurement weak. On the 
other end a measurement of the mixing fre- 
quency ratio ArUs/ Am^ would fix the ratio of 
hadronic parameters thus allowing a much tighter 
constraint on the unitarity triangle. The ex- 
pected range for Ar^s from the analysis of all 
available B and Kaon physics data within the 
SM is getting narrower and narrower (Arrig G 
[15.6,23.1] @ 95%C.L. [T]). Observing B° mix- 
ing outside the SM favored range would be a sig- 
nal for new physics. Experiments at LEP and 
SLD as well CDF in Run I have only provided 
a lower limit Am^ > 14.5@95% C.L. [2|. The 
high frequency expected for B° flavor oscillations 
poses formidable experimental challenges related 
to the need for extremely accurate proper time 
resolution in order to resolve the oscillation pat- 
tern. This is clearly reflected in the approximate 
expression for the significance of the mixing am- 
plitude measurement 3; 
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where S is the signal sample size, B is the back- 
ground, e is the efficiency for flavor tagging and D 



is the tagging dilution defined as = 1 — 

2w, where the mistag rate w — NiYf{Nji + Nw) 
is the fraction of B candidates with an incorrectly 
assigned flavor. The effect of finite resolution on 
proper time (ctct) reduces the statistical signifi- 
cance of the measurement, especially at high val- 
ues of Anris. A key point of the CDF upgrade for 
Run II has been the capability of triggering on 
hadronic decays of B mesons allowing the collec- 
tion of fully reconstructed B^ decays. Despite the 
limited statistics available in these modes the ex- 
cellent proper time resolution significantly helps 
the current sensitivity of CDF and will offer the 
unique chance of a determination of Arris from 
CDF data alone. We describe here the first search 
for _Bs mixing using 360 pb~^ of Run II data. 

2. CDF II Detector and Trigger 

The upgraded CDF detector for Run II 4 has 
been designed to enhance the sensitivity to 5° 
mixing, in particular with the LOO 5 layer of sil- 
icon detector very close to the beam pipe, the 
TOF system |S] and the Silicon Vertex Tracker 
(SVT) [Jj. This latter device, reconstructing on- 
line tracks in the silicon tracker, enabled trigger- 
ing on _B-decay vertices, for the first time at a 
hadron collider. The current CDF 6— physics re- 
lated trigger menu is based on the following se- 
lections: 

• Di-muons with muon transverse momenta 
Pt > 1.5 GcV/c to select b J/ipX decays. 

• Single muon or electron with Pt > 4 as- 
sociated with a displaced trigger track with 
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Figure 1. B°^D+TT-,D+~^(j)TT+ signal. 
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Figure 2. B'^.-^D+rX ,D+^(j)n+ signal. 



impact parameter greater than 120 fim (lep- 
ton+SVT). 

• A Two Track Trigger (TTT) requiring two 
displaced SVT tracks with impact param- 
eters greater than 120 fim, consistent with 
originating at a vertex 200 /im away from 
the beam axis. 

These different triggers provide large samples of, 
respectively, Bd,u^ J/'ipK^*\ Bd,u^lD'^*^X and 
Bd.u^ D'^*'^tt used to calibrate the vertex resolu- 
tion function, the flavor-tagging dilution and to 
extract the impact parameter trigger bias. These 
measurements are an essential part of the Bg mix- 
ing study making the latter really the "ultimate" 
6-physics analysis at CDF. 

3. Signal samples 

Currently we reconstruct B^ Dfn^ 
(hadronic analysis) and B^ D'^l~ X channels 
(semileptonic analysis). In both cases the i^^— > 
^7r+,£)+ -> K*°K+ and D+ 7r+7r-7r+ decay 
modes are used, where the light meson resonances 
are detected in the cji^K^K' and K*°^K+tt~ 
channels only. These modes combined correspond 
to 4.9% of the total Ds decay rate. 



i?^ meson signals are identified as sharp peaks 
either in the reconstructed B^ mass for hadronic 
decays or in the Ds mass for opposite sign com- 
binations of a lepton (/i or e) and a D^ meson 
candidate. 

B meson candidates are reconstructed starting 
from charged tracks and leptons matched to trig- 
ger objects satisfying on-line trigger selections, in 
order to accurately model trigger biases to the 
signal proper time distribution through Monte 
Carlo (MC) simulations. To ensure optimal ver- 
texing and momentum resolution only tracks re- 
constructed both in the central drift chamber and 
in the silicon vertex detector with a minimum 
transverse momentum of 400 MeV/c are used. A 
secondary/tertiary vertex fit is performed on all 
possible track combinations using a mass con- 
straint on Da mesons involved in hadronic decays. 
Further selection criteria have been designed to 
optimize S/ ^/{S + B) , using MC simulations for 
signal and sideband data for background. 

Hadronic decays are selected requiring a signif- 
icantly displaced secondary vertex and a tertiary 
vertex consistent with a downstream decay of a 
charmed meson. Selecting high probability ver- 
tex fit reduces combinatorial background and a 
minimum Pt on the pion from B^ decay is re- 
quired to reduce background from Ds production 
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Table 1 

B° yield and S/B in 360 pb^i of CDF II data. 



Decay Channel 


Yield 


S/B 






526 ± 33 


1.8 






254 ±21 


1.7 




7r~ , (tt^ TT^ 7r+ ) 


116 ± 18 


1.0 




■^D+l-X,i(j)Tr+) 


4355 ± 94 


3.1 


B^- 


-*D+l-X,iK*°K+) 


1750 ±83 


0.4 


B's- 




1573 ± 88 


0.3 



in pp—^ cX events combined with a fragmentation 
track. Remaining background sources are com- 
binatorics and partially reconstructed B decays. 
The former is modeled by exponential functions 
and fixed by fits to the high mass sideband 
while the latter appear as a well separated bump 
to the low mass side of the signal peak. From 
MC simulation we derive shapes for the contribu- 
tions from decays like 5°-^ L'*7r^ 1)5(7, 7r°)7r+ 
or i3g — > > Ds7r"'"(7r°) where the 7 or 7r° es- 

cape undetected. Other backgrounds considered 
are Cabibbo-suppressed B^-^DgK decays (whose 
contribution is fixed to the analogous i?" decay 
rate) and reflections from b-baryon and de- 
cays whose shapes and relative contributions are 
fixed to those expected from simulation. 

Similar selection criteria are used to isolate 
semileptonic decays. Further requests on the K* 
and (j) helicity angle and on the the lepton + Ds 
invariant mass (2.3 < mw^ < 5 GeV/c^) reduce 
combinatorial and cc background. An important 
background (of order 20% of the signal sample) 
arises from B^d\^^D^*)^X with D^*J^^lvX de- 
cays. The expected contribution is determined 
from MC using world average branching ratios 
and considered in the mixing analysis with proper 
lifetime and mixing probabilities. The back- 
ground from real prompt Dl ' paired with a ran- 
dom track faking a lepton (prompt background) 
is estimated from wrong sign, l^D^, signal. Its 
contribution is about 4% of the right sign sample. 

An example of the clean signals observed in 
hadronic and semileptonic decays is shown in 
Figs, n and [3 Number of events and S/B val- 



ues extracted from fits to the B^ or Ds invariant 
mass are displayed in Tabled 

4. Proper Time Reconstruction and Reso- 
lution 

The proper decay time of B'^ mesons is recon- 
structed using the distance in the transverse plane 
from the beam axis to the decay vertex (L^y) 
and the transverse momentum of the candidates: 
■ mBjPr- For semileptonic decays 
a K factor is needed to compensate for missing 
momentum carried by neutrinos ctb^ = K ■ ct^ 
where CTg is reconstructed from the lepton + 
Ds transverse momentum. The K factor distri- 
bution is derived from MC simulation of inclusive 
semileptonic decays. 

The time dependent mixing probability is con- 
voluted with a Gaussian resolution function us- 
ing the event-by-event expected error on ctb^ ■ 
The latter includes a scale factor derived study- 
ing the proper time resolution function on a sam- 
ple of pseudo-B^ candidates constructed from a 
prompt Ds meson paired with a fragmentation 
track. This sample exhibits characteristics very 
similar to those of the signal and allows the pa- 
rameterization of the scale factor as a function of 
several kinematic variables. The average ct un- 
certainty for the golden B^ ^ Dfir^ ,Df (j)TT^ 
mode is « 100 fs corresponding to 60-70 ^m res- 
olution on the Bs decay length. Note that the 
exponential term in eq.njiniphes a 60% loss of sta- 
tistical power for 100 fs resolution at Am^ = 15. 
Improving the proper time measurement accu- 
racy, e.g. by using an event-by-event primary 
vertex, will certainly help in a future updated 
analysis. 

5. Lifetime Measurement 

Determinations of i?-meson lifetimes in sev- 
eral exclusive B^ and B^ decay channels provide 
an important check on the accuracy of proper 
time measurement. This has been successfully 
performed [HI using B^ J/ipK'-*'' and semilep- 
tonic decays collected with a version of the trig- 
ger with no impact parameter requirement on 
tracks. For the present mixing analysis an ad- 
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Tagging dilution, efficiency and calibrated sD^ (%) with statistical and systematic uncertainty. 


Tagger 


V< > 


s 


eiDSofisemil.) 


e{DSDf{had.) 


SMT 


35.9 


5.0 


0.557 ±0.047 ±0.034 


0.46 ±0.11 ±0.03 


SET 


29.4 


3.5 


0.293 ± 0.033 ±0.017 


0.18 ±0.06 ±0.02 


JVX 


16.9 


9.8 


0.263 ±0.035 ±0.021 


0.14 ±0.07 ±0.01 


JJP 


11.5 


14.0 


0.150 ±0.026 ±0.015 


0.11 ±0.06 ±0.01 


JPT 


5.0 


52.0 


0.157 ±0.027 ±0.025 


0.24 ±0.09 ±0.01 


Total 


35.9 


5.0 


1.429 ±0.093 


1.12 ±0.18 ±0.04 



ditional crucial point has been the removal of the 
bias introduced by the trigger and offline selec- 
tion criteria. A detailed Monte Carlo simulation 
of the trigger, including the effects of the varia- 
tions with time of the beam axis position with 
respect to the detector and of trigger algorithms, 
has been used to derive the efficiency clS Qti func- 
tion of proper time for each of the decay modes 
analyzed. This method has been validated mea- 
suring _B+ and lifetime in Dn modes collected 
by the same trigger used for the mixing analysis. 
The results (t(B+)(ps) = 1.661 ± 0.025 ± 0.013, 
t(BO)(ps) = 1.511 ± 0.023 ± 0.013, t(BO)(ps) = 
1.598 ± 0.097 ± 0.017) are consistent with world 
averages. The systematics are at the level of 1% 
and demonstrate that the the trigger bias is well 
under control. A similar analysis is currently un- 
der way using the same semileptonic sample 
used here. When completed, it will give the best 
measurement of lifetime in a flavor specific de- 
cay. For the time being a result consistent with 
the world average with a statistical uncertainty 
of 2% has been reported by CDF. 

6. Flavor Tagging 

Identifying the initial flavor of the B mesons at 
hadron colliders has always been a critical task 
due to the busy environment of b events in pp 
collisions. To determine the B^ flavor algorithms 
have been developed relying either on the other 6- 
hadron in the event. Opposite Side (OS) taggers, 
or on the expected charge correlation [B^ — K'^) 
between the B^ flavor and the kaons emerging 
close to the B^ flight direction from the frag- 
mentation and hadronization process. Same Side 



Kaon (SSK) tagger. Five different OS taggers are 
used here: 

• A Soft Muon Tagger (SMT), based on the 
detection of muons, with Pt > 1.5 GeV/c 
and pseudorapidity \ri\ < 1.5 from b~^ii~X 
decays. The dilution is parameterized as 
a function of a muon likelihood and muon 
momentum relative to the jet axis, P^''', to 
account for the higher dilution expected at 
higher P^*' due to lower 6— >c— >Z+ sequen- 
tial background. 

• A similar algorithm, Soft Electron Tagger 
(SET) , uses b decays to electrons with Pt > 
2.0 GeV/c detected in the central electro- 
magnetic calorimeter. 

• Three different Jet Charge Taggers, where 
a weighted sum of track charges in oppo- 
site side jets is correlated with the initial 
6-quark charge sign. The three algorithms 
differ on the choice of the tagging jet: jet 
containing an identified secondary vertex 
(JVX), jet characterized by at least one dis- 
placed track (JJP) or the highest Pt jet in 
the event (JPT). 

All the taggers have been tuned on an inclusive 
and almost pure sample of « 10^ semileptonic 
b decays from the lepton + SVT trigger. The 
dilution of each different algorithm is determined, 
as a function of several variables, from the the 
trigger lepton charge after correction for mixing 
and b^c^l^ backgromid. Taggers are ranked 
in order of their expected dilutions and only the 
decision from the highest dilution tagger available 
in any single event is used. 
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Since the calibration sample for the dilution 
may be different in certain aspects (e.g. Pt 
spectra) from the signal sample a scale factor 
(Sd) is allowed for each of the 5 tagging algo- 
rithms. The Sd factors for the semileptonic and 
hadronic analysis are fixed by measuring the i?° 
oscillation frequency and its amplitude using re- 
spectively the high statistics Bd.u — > IDvX and 
Bd.u-^ J/ipK'^*'' /Bd^u-^ D^*''Tr samples. A sum- 
mary of the tagging figure of merit eD^ is re- 
ported in Tab. 121 The scale factors are consistent 
with 1 within the precision allowed by the present 
statistics. 

7. Results and Perspectives 

Following 3 an amplitude scan is performed, 
where, at each value of the mixing frequency, an 
amplitude A is fitted from the data. The am- 
plitude A should be consistent with for mixing 
frequencies far from the true value of Am^ and 
with 1 in its vicinity. This method allows for a 
straightforward way to combine several measure- 
ment (semileptonic and hadronic) as well as re- 
sults from different experiments. 

We determined systematic uncertainties on the 
amplitude A using large toy MC experiments to 
minimize the effect of the limited statistics. The 
biggest effects on systematic uncertainties at high 
value of Arrig come from the level and the dilu- 
tion parameterization of the prompt background 
for the semileptonic analysis (contributing aA — 
0.15 @ Arus — 20ps^^) and the dilution scale fac- 
tor for the hadronic analysis {aA — 0.3 @ Anis = 
20ps~^). It should be noticed that the latter will 
certainly be reduced with more data since it is 
limited by the available number of Bd^u hadronic 
decays. In any case the combined systematic un- 
certainties are much smaller than the statistical 
uncertainties on the present measurements, for 
any value of the mixing frequency. 

The present CDF sensitivity is dominated by 
the semileptonic result: Anis > 7.7 ps~^@ 95% 
C.L. with a sensitivity of 7.4 ps"^. The com- 
bined semileptonic and hadronic analysis limit is 
Arris > 7.9 ps~^@ 95% C.L. while the combined 
sensitivity increases to 8.4ps~'^. 

Combining the CDF mesaurement with previ- 
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Figure 3. Amplitude scan for the combined B'^- 
D+l-X and B^-^D+ir- analysis. 



ously available results does not change the actual 
limit and only marginally improves the sensitiv- 
ity (from 18.2 to 18.6). It can be noticed, how- 
ever, that the current sensitivity of CDF II data 
is relatively better behaved at higher Am^ than 
most of previous experiments due to the valuable 
hadronic sample that will dominate the sensitiv- 
ity in future higher statistics searches for os- 
cillations. 
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